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We assess the effects on Big Bang Nucleosynthesis (BBN) of lepton number generation in the 
early universe resulting from the two-doublet four-neutrino mass/mixing scheme. It has been ar- 
gued that this neutrino mass/mixing arrangement gives the most viable fit to the existing data. 
We study full 4x4 mixing matrices and show how possible symmetries in these can affect the 
BBN 4 He abundance yields. Though there is as yet no consensus on the reliability of BBN calcu- 
lations with neutrino flavor mixing, we show that, in the case where the sign of the lepton number 
asymmetry is unpredictable, BBN considerations may pick out specific relationships between mixing 
angles. In particular, reconciling the observed light element abundances with a ^ v £ oscillation 
interpretation of LSND would allow unique new constraints on the neutrino mixing angles in this 
model. 

PACS numbers: 14.60.Pq; 14.60.St; 26.35.+c; 95.30.-k 



I. INTRODUCTION 

In this paper we revisit the long standing issue of how neutrino flavor mixing in the early universe might affect Big 
Bang Nucleosynthesis (BBN). Though this is an old problem |l],|2|,|3|,[I|,[| , with the advent of modern experiments {e.g., 
solar, atmospheric, and accelerator-based oscillation experiments) we can hope to begin constructing the neutrino 
mass/mixing matrix. Here we adopt the leading model for this, the two-doublet four- neutrino mass scheme We 
discuss this mixing in detail, going beyond recent good assessments of this problem [|]|],[l(J , to examine what BBN 
considerations ultimately may be able to tell us about parameters in this mixing scheme. 

Historical attempts []l],^]|J],[j,|5| at constraints on neutrino mixing are based on the argument that an active-sterile 
neutrino mixing that is too large at (or prior to) the BBN epoch will populate the sterile neutrino sea. The resultant 
increase in the total energy density of the universe at a given temperature speeds up the Hubble expansion. In turn 
this leads to a higher weak-freeze-out temperature and, consequently, it could lead to a higher neutron-to-proton ratio 
at Nuclear Statistical Equilibrium freeze-out (however, the neutron-to-proton ratio is determined not only by 
the expansion rate, see below). Since essentially all neutrons are incorporated into He in the early universe, such a 
mixing yields a higher 4 He abundance Y. This potentially could contradict the observationally inferred primordial 
4 He abundance. From the inferred abundances of 4 He and D/H there are strong constraints on the increase of the 
predicted Y from an increased energy density due to sterile neutrino production. These constraints are often translated 
into a limit on the "effective number of neutrinos" of 

1.7 < Nf < 3.2 (1) 

at 95% confidence that is strictly limited by the baryon-to-photon ratio, ?y, determined by the inferred relative 
abundance of deuterium (D/H) [Q. This limit is shown in Fig. [l], where we have calculated Y and D/H over a range 
in and i]. We have included corrections to Y derived from the zero and finite-temperature radiative, Coulomb and 
finitc-nucleon-mass corrections to the weak rates; order-a quantum-electrodynamic correction to the plasma density, 
electron mass, and neutrino temperature; incomplete neutrino decoupling; and numerical time-step effects p3[ . The 
observationally-inferred primordial mass fraction we adopt is generous: 0.228 < Y < 0.248 (95% confidence range) 
Q; and the inferred relative abundance D/H we adopt is the well-established D/H w (3.4 ± 0.5) x 10~ 5 |Q. For 
arguments and observational evidence on the reliability of the deuterium abundance see Ref. |Tq). Similar limits on 
N^ s have been derived by other groups fll6| , p"7 

neutrino sea (iV° ff = 4) in the BBN epoch. The concept of N^ n (an expansion rate measure) as the sole determinant 
of the 4 He yield is a misleading and dangerous one. 

A few 4x4 models have been advanced as simultaneous explanations of the atmospheric neutrino data, the solar 
neutrino data, and the data from LSND ^fy. One hierarchical scheme would have (near) maximal ^ v r mixing 
for the atmospheric deficit, a v e ^ mass splitting that gives either a vacuum or MSW solar solution, and the 
LSND indication of a large Sm 2 is given by "indirect" mixing through a sterile. This is Scheme I in Fig. 0. In this 



18 1 . These limits preclude by a wide margin a fully populated sterile 
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scheme, the mass states most closely associated with the active neutrinos form a triplet that has a significant mass 
splitting from the mostly sterile mass state. This scheme is ruled out by BBN because it requires large # v s and 
v e ^ v s mixing amplitudes to explain the LSND results through the indirect conversion — > v s — > v e (and the 
CP conjugate). The required large mixing amplitudes for Vj, ^ v s and v e ^ v s populate the sterile sea in the early 
universe through direct oscillationproduction of steriles P,^,BM] . This mass model is also disfavored by a combined 
4x4 experimental data analysis jjj ■ 

An alternative two-doublet hierarchical scheme (Scheme II) has (near) maximal ^ v s mixing to fit the atmo- 
spheric data in Super-Kamiokande, and v e ^ v T mixing explaining the solar neutrino puzzle. This option is excluded 
by BBN since the ^ v a transformation is too large to be compatible with 4 He observations |^],[l9| ■ (Recently Foot 
and Volkas [p3[ have argued that ^ v s maximal vacuum mixing for SuperK atmospheric neutrinos in a related 
mass scheme can be reconciled with BBN limits.) El 

Therefore, adopting the previously considered limits from BBN, we are left only with the two-doublet four-neutrino 
mixing model, Scheme III. The BBN effects of the model in Scheme III were considered in Ref. ||. In this paper we 
expand on the analysis of the mixing and suggest how BBN could give potentially stricter limits that those found in 
Ref. d). (We note here that the "inverted scheme," where the v e /v s doublet is more massive than the v^/vr doublet 
is not yet completely ruled out by laboratory and astrophysical considerations.) 

A new twist was added to the saga of neutrino-mixing in the early universe when it was found that resonant active- 
sterile neutrino transformation in the BBN epoch can alter neutrino energy spectra and generate a lepton number 
asymmetry f2"l"| , |22|| . This raises the possibility that an asymmetry in v e /v e numbers could be generated. In the mean 
time, their energy spectra may be modified so that the weak reaction rates themselves may change, resulting in a 
different neutron-to-proton ratio and a different 4 He yield Y |23|,^4|,^5|] . It has been argued that a positive v e jv e 
asymmetry behaves like a positive chemical potential for v e . This would reduce Y, and by the same token a negative 
asymmetry would increase it [ p3| . However, this argument of a direct asymmetry- V leverage relation is too naive in 
the context of active-sterile neutrino mixing, and is in fact incorrect. This is because the process of lepton number 
generation via resonant active-sterile neutrino mixing potentially has a crucial and unique feature [ p7| , p2|j2^j27| , |28t : 
that the lepton number asymmetry is first damped to essentially zero and then can oscillate chaotically with an 
increasingly larger amplitude, until it converges to a growing asymptotic value that is either positive or negative. 

However, the existence of "chaoticity" (unpredictability) in the final sign of the lepton asymmetry, L, is controversial. 
Indeed, all claims (e.g. pl"| , p2[ ) for lepton number generation in the early universe via neutrino mixing are now at 
odds with at least one calculation (Dolgov, Hansen, Pastor and Semikoz [p9|| ). Ref. [[29j uses a new formulation of 
the solution of the neutrino energy density matrix evolution that finds a non-chaotic generation of only a small and 
insignificant lepton asymmetry. Refs. |2(| and (27) recalculated the lepton number generation found in Ref. |^2| and 
corroborated a random nature to the sign of the L. Ref. ]27j ] showed that certain over-simplifying approximations 
in Ref. |29| may have unphysically stabilized the evolution of the lepton number. Ref. [p8[ found a randomness in 
the sign of L, but only for small non-phenomcnological mixing angles, much smaller than those that may be involved 
in neutrino oscillation solutions to the experiments discussed here. The oscillations in lepton number sign seen in 
Ref. |28|| only occur below the precision of their numerical solution. One may call into question then the ability of 
their numerical formulation to resolve oscillations occuring at larger mixing angles which can be found in the more 
straight-forward momentum-averaged solution of Refs. |22| , |26| , ^7|] . The momentum-averaged solution is valid at the 
relevant high-temperatures where the instantaneous approximation to repopulation are good. For the LSND mass 
scale in Scheme III considered here, m 2 — m 2 e w to 2 — m 2 s ~ 0.2 to 10 eV 2 and a small effective mixing angle 

(sin 2 8 c s > 10 -10 ), the bulk of the active neutrinos undergo resonance at a temperature T ~ 10 to 20 MeV where the 
instantaneous approximation can be considered valid. 

Given the controversy and disagreement among these different calculations it is difficult to take any BBN-derived 
constraint with confidence. With this caveat in mind, we explore what constraints might be possible if the lepton 
number generation magnitude is as in Refs. [^1 22 1, but where chaoticity in lepton number sign obtains for relevant 
parameters. 

In a chaotic lepton number generation regime, the sign of the lepton number asymmetry is independent of the 
initial conditions before the amplification begins, and is exponentially sensitive to the parameters involved during the 
chaotic oscillitory phase |22j . In turn, the sign of the asymmetry cannot correlate over a scale bigger than the particle 
horizon jnj]. 



1 There is still another loop-hole: if there is a pre-existing lepton number asymmetry with a magnitude > 10 5 during the 
BBN epoch, the ^ v s mixing can be suppressed This assumption must then involve physics that is beyond neutrino 
mixing. 
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This causal structure of space-time can make it impossible to obtain a universe with a uniform lepton number 
asymmetry. Instead, the lepton number generating process gives rise to a universe with numerous lepton number 
domains with similar lepton number magnitudes, albeit different signs p9| . The size of each domain is less than the 
horizon size at that epoch (~ 10 10 cm at the weak-freeze-out temperature, although the detailed geometric structure 
of the domains at the BBN epoch depends on the manner in which domain "percolation" occurs) . The distribution 
of domains with different signs is completely random so that in total each sign occupies half of the space. 

The overall primordial 4 He yield in such a universe must be an average of Y over domains with opposite v e /v e 
asymmetries. Interestingly, this implies that the overall Y is always larger than that expected when there are no 
lepton asymmetries. This is because the increase in Y from a negative lepton asymmetry generated by the resonant 
active-sterile neutrino mixing process is always bigger than the decrease in Y from a positive asymmetry generated 
in the same process Q . 

In this paper we examine in detail how this causality effect operates in the context of a specific four neutrino mixing 
scheme. There are some surprises. 

We can quantify these arguments for the two-doublet neutrino mixing model (Scheme III). In this model, the 
Vn v e mixing that fits the LSND data has a mass-squared-difference — m 2 c ~ 0.2 to 10 eV 2 and an effective 
two-neutrino mixing angle satisfying sin 2 20 eM ~ 10~ 3 to 10 -2 |5(J; the v e ^ v s mixing that solves the solar neutrino 
puzzle has either m 2 — m 2 ~ 10 -5 eV 2 and an effective two-neutrino mixing angle satisfying sin 2 28 es ~ 10~ 2 
for the small mixing angle (SMA) MSW (Mikheyev-Smirnov-Wolfenstein) solution or |m 2 — m 2 | ~ 10~ 10 eV 2 and 
sin 2 26 es ~ 1 for the vacuum solution J3l| ; the Super-Kamiokande atmospheric neutrino data are nicely fit by ^ v T 
mixing with |m 2 — to 2 t | ~ 10~ 3 to 10 -2 eV 2 and sin 2 26 ^ T ~ 1 |p2| . The hierarchy in masses which emerges from 
this mixing scheme includes an upper doublet of heavier neutrinos, consisting of (almost) degenerate Vu and u T , and a 
lower doublet of lighter neutrinos, consisting of slightly mixed (in the case of the MSW solution) or (near) maximally 
mixed (in the case of vacuum mixing) u e and v s . The inter-doublet mixing between and v e is small. So also must 
be the inter-doublet mixing between or v T and i/ s , so as to avoid conflicts with BBN ^^g- 

The large mixing angle (MSW) solution to the solar neutrino problem has been previously disfavored by BBN |3|,|],|| . 
We consider the SMA MSW active-sterile neutrino mixing solar solution and the vacuum active-sterile neutrino mixing 
solar solution. In Ref. pi]] , the vacuum sterile neutrino solar solution was considered to be disfavored using the current 
combined solar neutrino experiment rate data. However, in Ref. |3j|, it was argued that if only the solar Super-K 
spectrum or the seasonal solar data are considered, then the vacuum active-sterile neutrino mixing solution gives a 
better fit than a vacuum active-active neutrino mixing solution. Since the nature of neutrino mixing of solar neutrinos 
is still not certain and since BBN considerations may prove to be enlightening, we entertain the possibility of a vacuum 
active-sterile neutrino mixing solution to the solar neutrino deficit. The vacuum sterile neutrino mixing solution to 
the solar neutrino problem was also considered in terms of a neutrino mass model represented by the [SU(3)] 3 or E 6 
groups in Ref. 34 1. 

The interpretation of current experimental results usually is framed in terms of an effective two neutrino mixing 
scenario (i.e., in terms mass-squared-differences and mixing angles) for each experimental situation. This model is 
approximately valid in the two-doublet four-neutrino mixing scheme because, as a result of the mass hierarchy, one 
two-species mixing dominates in each of the above experiments. It is, however, more informative to employ the full 
4x4 mixing matrix in our discussion. In the next section we will briefly review what has been learned about this 
mixing matrix from the current experiments. We will then proceed to consider BBN 4 He synthesis in the presence of 
hierarchical four-neutrino mixing schemes. From the 4 He yield we infer potential new constraints on the inter-doublet 
mixing matrix elements between active neutrinos and the sterile neutrinos. In section III, we will summarize our 
results. 



II. THE HIERARCHICAL FOUR-NEUTRINO SCHEME AND THE PRIMORDIAL 4 HE ABUNDANCE 

We adopt the convention of employing Greek indices to denote flavor eigenstates v Sl v e , and v T , and employing 
Latin indices to denote mass eigenstates z/o, fx, V2 and ^3. The two bases are related by a unitary transformation U : 

V a = ^ U ai v i- ( 2 ) 

z=0,3 

The mass matrix in the flavor basis is then 

M af} = E E K k m k S kl U^, (3) 

fc=0,3 (=0,3 
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where rrik are the mass eigenvalues, and 5ki are the Kronecker deltas. In the scheme considered here, mo, mi <C 
m 2 , m 3 . 

The full expression for U can be found, for example, in Eq. (8) of Ref. |3jJ . It has 12 degrees of freedom, parametrized 
by 6 rotation angles 9 a b, and 6 CP- violating phases 4> a b- The parameters are symmetric in indices a and b, which 
run from to 3 (a < b). We follow convention and write s a b = sin^e"^ 6 and c a b = cos 9^. Because the inter- 
doublet mixing is small, \U s2 \ 2 , \U s3 \ 2 , \U e2 \ 2 , \U e3 \ 2 , |[/ p o| 2 , |l^i| 2 > \ u ro\ 2 and \U T i\ 2 are small, implying that 
| S02 1 , | S03 1 , I si2 1 , I S13 1 are small. In fact, the Bugey result can be translated into a limit |si2|, | S13 1 < 0.1 |J5). The 
assumption that and v T are (nearly) maximally mixed yields C23 ~ | S23 i ~ l/y/2. Furthermore, the LSND result 
suggests that it is likely that | S12 1 , | Sis) ~ 0.1 [^5|. Finally and obviously, we should have C02 ~ Q53 ~ en ~ C13 ~ 1. 

To leading order in S02, so3j S12, and S13, we have |35| 

C01 SqI s 02 s 03 \ 
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-S23 C23 



(4) 



/ 



In Ref. |35| this matrix is discussed and is shown to be unitary to second order in the LSND mixing angle. 

In a simplifying approximation, one can take all mixing angles except #01 > ^12 and #23 to be zero, ignore the CP 
violating phases, and take the Super-K associated mixing to be maximal, C23 — l/y/2. In this case, the transformation 
matrix becomes 



U 



( 



coi 
-soi 
S01S12/V2 
-S01S12/V2 



Soi 
Coi 

-C01S12/V2 
C01S12/V2 





S12 

1/V2 
-1/V2 



which is (approximately) unitary for nearly maximal mixing, C23 — 1/V2- 
We can define a linear row transformation P6Lfi7| 



and 



V2 




(5) 



(6) 



such that Eq. (0) becomes 
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(9) 



Here, we can see that the fourth state \v*) is a mass eigenstate. If there are no lepton asymmetries generated by a 
mechanism other than neutrino mixing, then the muon- and tau-neutrino flavors see the same matter effects (that is 
the same thermal and fermion potentials) throughout their evolution. The state \v*) is mixed with the sterile and 
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electron neutrino, and will undergo resonant MSW transformation under the appropriate conditions. However, \v*) 
will pass through resonances unchanged. This reduces the 4x4 mixing to essentially a 3 x 3 evolution, at least as far 
as MSW resonances are concerned^ It should be noted that the mixing matrix discussed here uses the rotation order 
convention U = R23R13R03R12R02R01, while Caldwell, Fuller and Qian use E/cfq = R23R01R12 [pql - The different 
unitary mixing matrices are physically equivalent under the same approximations, and exhibit the same decoupling 
seen in Eq. (||). 

The matter effects present at the epoch of BBN comprise two pieces: a part due to a finite-temperature thermal 
bath, and a part due to a possible lepton number asymmetry in the active neutrino sectors (the baryon number 
asymmetry and the associated electron-positron asymmetry is too small to play a significant role in the neutrino 
mixing). For two-neutrino mixing, the effective matter mass-squared-difference is 



5m 2 a (3 (eS ) = \ (ml a - ml Y sin 2 26 af 



1/2 

(10) 



and the effective mixing angle satisfies 

fl (cff.) ( m L -m^)sin2g aA 

tanf ^ ~ (ml a -ml fi )cos29 a , + 2EV^ + 2EV a y ' ' 



Here we use m Va to denote the mass eigenstate most closely associated with a neutrino of flavor a; E is the neutrino 
energy; is the effective potential due to the finite-temperature part of the matter effect, and V^a is the effective 
potential due to the lepton number asymmetry. Both potentials vary with temperature T. 

We will first consider the case that the lepton number asymmetry is negligible, 2EV£g <C 2EV^, (m 2 a — 
m 2 ) cos 29 a p. Mixings between active neutrino species alone do not modify the 4 He synthesis because active neutrinos 
share the same number density distribution in the BBN epoch. (This is not rigorously true because electron/positron 
annihilation overpopulates v e jv e slightly but its impact on the 4 He yield is less than 0.1% PH .) Mixings between 
active neutrinos and sterile neutrinos, however, convert active neutrinos and so populate initially unoccupied sterile 
neutrino states. These mixings can therefore affect the energy spectra of active neutrinos, as well as the 4 He yield. In 
particular, MSW resonances, at which the local effective mixing reaches maximal values, can occur between an active 
neutrino species v a and the sterile neutrino species when (rn 2 ,^ — m 2 v ^) cos2# QS + 2EV^ S — 0. Since 

V T „ _ A n Va +n Da ^ ^ 
n 7 

where A w 105(30) for a = e (fi, t), this resonance condition is met when the temperature of the universe is 

-1/3 I (rm 2 _2 ^ „„„oa I/ 6 



To 



(m 2 — m 2 , ) cos 29 D 



leV 2 



(13) 



where T » 19(22) MeV for a = e (n, r). 

Therefore, for a two-family mixing between v* and v s , with m 2 ,. — m 2 s ~ m 2 , — m 2 ^ ~ 0.2 to 10 eV 2 and a 

small mixing angle, the bulk of the active neutrinos undergo resonance at a temperature T ~ 10 to 20 MeV. This is 
long before v* and v* decouple thermally/chemically from the thermal background. As discussed briefly earlier, the 
v* ^ v s mixing potentials behave identically to the v* v s in the BBN epoch, and they undergo the same scattering 
and collision evolution. Thus, the approximations giving Eq. (Q) effectively "decouple" v*, with only u* going through 
resonances. In the following discussion, v* decoupling can be assumed; however, even if the approximations leading 
to Eq. (H) are invalid so that angles other than 9q\^ 9\2 and #23 are non-zero, the following discussion is still relevant 
for the standard and can be extended to v T . 

Resonances may also occur for v e v s mixing with m 2 — m 2 ~ 10~ 10 eV 2 (a possible vacuum solution for the 
solar neutrino problem). But the resonance temperature is < 0.01 MeV, which corresponds to an epoch long after 
the weak-decoupling of neutrinos. Resonant active-sterile neutrino mixing within the lower doublet in this case then 
cannot influence BBN. Therefore, the BBN constraints on mixing models involving a vacuum solar solution are less 
stringent. 



2 Note that Ref. Ji| also points out that one "linear combination" of and v T "oscillates" with v a while the other decouples. 
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In the two-doublet neutrino mixing schemes considered here, the v* v s channel is essentially decoupled (at or 
near its resonance temperature) from the other mixings within the four species family. Therefore, we can take the 
two-neutrino mixing picture as applicable. This is so because while the v* # v s channel is matter-enhanced at its 
resonance, the other mixings are not, or even are suppressed by the matter effects by a factor of |(to 2 q — ml^) /2EV^\ 2 
with respect to their vacuum mixing amplitude. For example, the inner-lower-doublet v e ^ v s mixing is suppressed 
by a factor ~ 10 9 for the SMA MSW mixing solar neutrino solution or ~ 10 19 for the vacuum mixing solar neutrino 
solution, because 2EV es ~ 8EV^ S w 4(m^ - m 2 J. Also, since 2EV fie = ^(V^ - V es ) w -SEV^ w 3(m£ -m£j « 
3(m 2 — to 2 b ), the inter-doublet ^* # ^ e mixing is suppressed by a factor ~ 10 with respect to its already small 
vaccuum mixing amplitude (< 10~ 2 ). Therefore, it is safe to employ the two-family mixing picture to investigate the 
v* ^ channel at or near its resonance. 

As pointed out in several previous papers, there are two possible consequences of a resonant active-sterile neutrino 
mixing: (1) the total neutrino energy density at a given temperature increases if neutrino pair production is still 
effective in replenishing the converted active neutrinos; (2) a lepton number asymmetry may be generated in the 
active neutrino sector from initial small and negligible statistical fluctuations during the resonant active-to-sterile 
neutrino conversion process. In either the nonresonant or resonant mixing case, the limit on the total neutrino energy 
density from the primordial 4 He abundance puts the following constraint on the parameters of active-sterile neutrino 
mixing |§: 

/ 2 2 \ • 4 a <r f 10~ 9 eV 2 if v a = v e ; , , 

{m Va - mj sin 26 as < j lQ _ 7 ^ . f ^ = ^ ^ (14) 

In the resonant case, the increase in the neutrino energy density is significant to BBN when the resonance is 
adiabatic and occurs before chemical decoupling of the active neutrino (about 5 MeV for i/*). The lepton number 
asymmetry generated from an initially very small asymmetry by the resonant active-sterile neutrino mixing process 
is significant to BBN if part of the asymmetry resides in the u e /P e sector and is of order > 0.01. In the two-doublet 
neutrino model, this is achieved by having a resonant v* ^ v s mixing generate L„» , and having a resonant v* ^ v e 
mixing transfer part of L„* into L u . 

In Ref. p4| , we calculated the change in the primordial 4 He abundance due to such a process. In the regime where 
10 _1 eV 2 < rn 2 * — m 2 < 10 eV 2 and sin 2 29 fl * s > 10~ 10 , the mixing angles are large enough to generate a lepton 
number asymmetry ^2^. Mixings with smaller mixing angles cannot have a material effect on BBN. In regions of the 
universe where L v * and in turn L„ e is positive, the n— >p rate is enhanced while the p— >n rate is reduced. This change 
in n^p rates tends to lower the neutron-to-proton ratio and consequently Y . 

On the other hand, the increased neutrino energy density from the active-sterile neutrino mixing before neutrino 
chemical decoupling always tends to increase Y, The overall result is a decrease in Y, as the former effect (v e /v e 
asymmetry) dominates. In places where L„» and L Vc are negative, however, the n^p rates are changed so as to 
increase Y, This is, of course, in addition to the increase in Y from the energy density effect. 

When averaged over the positive domains and the negative domains, the net Y turns out to be consistently larger 
than that predicted by the standard BBN picture assuming no neutrino mixing. This is because the increase in Y in 
negative domains is never compensated by the decrease in Y in positive domains, as shown in Fig. |^ |22] , [l^] . 

This result is rather different from that in Ref. || and given the controversy in the calculations alluded to above, it is 
difficult to say which, if either, is correct. Part of the discrepancy may be because the effect of energy density increases 
on Y were not appropriately taken into account in Ref. || . A simple average of the change in Y in our calculation 
is not beyond observational uncertainty bounds. The average AY" is always less than ~ 0.001, or, cquivilantly, N^f 
is always less than 3.08. If a definitive, confident solution to lepton number generation by neutrino mixing in the 
early universe were to show unambiguously that the sign of the neutrino asymmetry for this specific range of neutrino 
mixing parameters is positive (negative) then the predicted change to Y would follow the lower (upper) curve. The 
positive lepton number result alone does not exceed the observational bounds of Fig. |l[ The observations more 
greatly constrain increases to Y. A negative lepton number result alone therefore creates a AY that is too large for 
<5to 2 , £ > 2.5 eV 2 . 

The averaged Y is only a lower limit to the actual Y in the two-doublet neutrino model with chaotic lepton number 
generation. This is because, as first pointed out by Shi and Fuller ]l9| , additional increases in Y may arise from 
an extra channel for increasing the total neutrino energy density. In the case where the sign of the lepton number 
generated is undetermined, the extra channel for sterile neutrino production results from the lepton number gradients 
between domains of regions with opposite lepton number sign (which are sub- horizon scale at the BBN epoch). The 
domain boundaries can meet the conditions for resonant conversion of not only v* to v s (and v* to v s ) but also v e to 
v s (and v e to v s ). Therefore, the sterile neutrino sea is not only populated within domains by the resonant neutrino 
mixing that drives the lepton number generation in the first place, but can also be populated at domain boundaries 
by the same resonant neutrino mixing as well as other active-sterile neutrino mixings. 



G 



To avoid fully thermalizing the sterile neutrinos, therefore, requires that this extra channel of sterile neutrino 
production be suppressed. In other words, all resonances of active-sterile neutrino conversion have to be non-adiabatic 
at domain boundaries. This yields another limit for the two-doublet neutrino scheme fjj| : 

sin 2 < 10~ 10 , for an MSW solution to the solar neutrino problem; 

(15) 



^ro 2 . — to 2 sin 2 28 p* s < 10 4 eV 2 , for a vacuum solar neutrino solution. 



Beyond these limits, Y increases by at least 0.013 due to a fully populated sterile neutrino sea resulting from active- 
sterile neutrino conversions at domain boundaries. 

These limits are summarized in Table 1. Because of the decoupling of various mixings, the above limits can be 
directly translated into constraints 

| s 2 1 , |so3 1 ^ 10~ 5 , for a SMA MSW solar neutrino solution 

or when m 2 « — ml > 4 eV 2 ; 



|*oa|, M < 10- 



-2 



(m 2 , - to 2 s ) / 



leV 2 



-i/2 (16) 
, for a vacuum solar neutrino solution 



and m 2 . - ml < 4 eV 2 



The constraint on the large mass difference is the result of the fact that for to 2 . — to 2 > 4 eV 2 domains of L„* 
will facilitate v* v s population of the sterile sea across domain boundaries. These constraints imply that the 
inter-doublet mixing elements of or v T with v s (sq2, S03) are ~ 10 4 times smaller than the elements associated with 
mixing with v e (S12, S13), if the solar neutrino problem has its roots in v e ^ v s mixing with to 2 — to 2 ~ 10 -5 eV 2 
(the SMA MSW solution). However, the mixings between the upper doublet neutrinos and the lower doublet v e and 
v s are within an order of magnitude if a vacuum "just so" v e # v s mixing with |to 2 — to 2 | ~ 10 -10 eV 2 explains the 
solar neutrino data. 

Let us now confine our attention to the latter scheme, where to 2 . — m 2 ^ is < 4 eV 2 and the solar neutrino deficit 

is explained by (nearly) maximally mixed v e # v s vacuum "just so" oscillation with 5to 2 s ~ 10 -10 eV 2 . This scheme 
can produce a significant increase in Y when the values of the inter-doublet mass splitting are in the range that is 
required to explain the LSND results. The primordial 4 He yield in this scenario is sensitive to the relative level of 
mixing for v^. T ^ v e and for v^ r ^ v s . For example, we can define a factor 

sin 2 20 Me 

— — 2 

sin 29^ s 

= \Sl2C23 + S13S23I 2 / 17 N 

|s 2C23 + S03S23I 2 ' 

where 0„ e and 6^ s are the effective two-neutrino vacuum mixing angles (to be constructed from the matrix elements 
in Eq. (Q)) corresponding to ^ v e and v s , respectively. Figure ^ then shows limits from the observed 4 He 

abundance for three cases: (1) F — 1; (2) F — 10; and (3) F — 100. More specifically, if the atmospheric neutrino 
problem solution is maximal ^ v T mixing (and non-CP violating), then F is just the ratio of inter-doublet mixing 
angles 

FJ-^±^ 2 . (18) 

|S02 + S03| 

It is noteable that a novel solution to r-process nucleosynthesis in Type II supernovae by Caldwell, Fuller and Qian 
fsij involves a 4 x 4 model where only S12 is non-vanishing in the above expression for F. In this case, effectively 
F — * 00, indicating an exceptionally high degree of "symmetry." (Here by "symmetry" we mean that only one mixing 
angle governs the inter-doublet mixing, not four; however, this is very asymmetric as far as F is concerned!) Both this 
solution to r-process nucleosynthesis and the domain-conversion-based BBN considerations discussed in this paper 
favor a large F. The mixing matrix in Ref. [^6| , which is related to that shown in Eq. |^, exhibits the symmetry 
that allows a decoupling of a neutrino state {\v*)) from MSW resonant evolution. For mass Scheme III to both 
account for the LSND signal (whose two- neutrino parameter space is shown in Fig. || |59|) and be consistent with 
domain-conversion-driven BBN effects, the difference in magnitude of the inter-doublet mixing angles must be large 
(a large value of F). 

If one were to entertain models where F was not large (i.e., comparable values for ^12,^13,^02 and #03), then 
already we can see from Fig. [| that the LSND data compatible with smaller F all lies near to 2 — to 2 w 1 eV 2 and 



7 



sin 2 2# Me w 10~ 3 . Concommitantly, future experiments such as BooNE |40) might indicate mixing parameters which 
fall outside this limit. Assuming our BBN domain-conversion-driven considerations are correct, this would argue 
strongly for large F. 

This would be a remarkable outcome. Naively, one might assume that, e.g., the BooNE experiment measures only 
the effective two-neutrino mixing 9^ e . However, when combined with other experiments and the BBN physics adopted 
here, it is evident that many other mixing matrix elements (those in F) are also probed. 



III. SUMMARY 



We have explored the possible change in the primordial 4 He abundance in the currently favored two-doublet four- 
neutrino {vT-v^jvz-Vs) mixing scheme proposed to simultaneously explain the current neutrino experiments. Though 
definitive calculations of matter-enhanced neutrino conversion effects arc elusive at present (in the eyes of some), here 
we have adopted the set of calculations which could give the tightest constraints on the neutrino mass and mixing 
matrix. We do this in the spirit of determining what may be possible. When we analyze the BBN effects in the 
context of full four- neutrino mixing we find some remarkable hints. Namely, we find that putative limits on the 
matrix elements governing v^. T ^ v e and f MjT ^ v s could be very restrictive. 

We have found that these limits strongly depend on the mixing between v e and v s : the limits are exceptionally 
strong if the v e # v a mixing parameters are in the range of the SMA MSW solution of the solar neutrino problem; 
but they are much less so, allowing the v^, v T ^ v a mixing to be at the same level as the v^, v T ^ v e mixing, if the 
v e v s mixing parameters lie at the vacuum "just so" solution region the solar neutrino problem. (Note: Since the 
submission of this paper there has been recent evidence from Super-Kamiokande that disfavors all sterile neutrino 
solutions to the solar neutrino problem jtoj . This emphasizes BBN constraints on any four- neutrino models.) In 
addition, we have found that if m 2 — rn~ « m 2 — m 2 > 4 eV 2 , the v^, v T x± v s mixing should also be extremely 
small. 

Therefore, the mixing structures in these different cases of the hierarchical four-neutrino schemes can be very 
disparate. Potentially unless the inter-doublet active-active and active-sterile mixings are very asymmetric, BBN 
considerations demand a vacuum "just so" solution to the solar neutrino problem and an LSND solution with m 2 — 

m 2 ~ 1 eV 2 and sin 2 28^ e w 1CP 3 in the two-doublet hierarchical mass scheme. Alternatively, future # v e 
experiments such as BooNE could be able to place significant constraints on functions of #12, #13, #02 and 6*03- This is 
a tantalizing result. Realizing it depends on the veracity of the particular BBN calculations we have adopted, but it 
is clear that the stakes are high. 

X.S., G.M.F. and K.A. acknowledge partial support from NSF grant PHY98-00980 at UCSD. K.A. wishes to 
acknowledge a NASA GSRP Fellowship. 
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TABLE I. BBN Limits on the inter-doublet active-sterile mixing in the two-doublet neutrino scheme. 





Type of Solar Neutrino Solution 


Limit on Effective sin 2 26 ^ s and sin 2 26 TS 


>4 eV 2 


SMA MSW 
Vacuum 


< lO" 10 


< 4 eV 2 


MSW 
Vacuum 


< 10" 4 eV 2 (mj„ ~miy L 




FIG. 1. The limits in the effective number of light neutrinos, N„ s , in BBN for varying barvon-to-photon ratio, rjio = i) 10 10 . 



The solid contours are 95% confidence limits on the inferred 4 He mass fraction, Y [)l4|, and the dashed contours are 95% 
confidence limits on the inferred relative abundance of D/H 
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Scheme I 



Scheme II 



Scheme III 



i/^i/g =: = : = : = : u^, u T = 

! ! ! 

^ TO LSND ^ TO LSND ^ m LSND 

I I I 



l/ e , V T 



FIG. 2. The three general mass heirarchies discussed. In all cases, the mass splittings correspond to 
<5"1lsnd ~ leV 2 ,<5m 2 tm ~ l(T 3 eV 2 , and <5m 2 olar ~ 10" 5 (10~ 10 ) eV 2 for the MSW (vacuum) solar solutions. Scheme I, 
its mirror (m„ c , m„ M , m„ T > m s ) and Scheme II (and its mirror) have previously been ruled out by BBN. In this paper, we 
consider constraints on mass heirarchies and mixings in Scheme III. 
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FIG. 3. The increase in the primordial 4 He yield in the two-doublet mass scheme, as a function of the inter-doublet 
mass-squared-difference. The mixing amplitude between v a and the v^-v T doublet is assumed to be not too small, > 10 -10 . 
The dashed curve: the increase in Y in individual domains. The solid curve: the increase in Y averaged over positive and 
negative lepton number domains. 
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sin 2 20 



FIG. 4. The effective ^ v e mixing parameters suggested by BBN and LSND for varying values of the asymmetry factor, 
F. Shown are the 90% and 95% CL limits for LSND. Curves labelled K, B, and N are the 90% CL limits from KARMEN2, 
Bugey, and NOMAD, respectively. Experimental confidence regions are adapted from Ref. |S9[. 
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